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ABSTRACT 

We report 60 and 90 fim observations of 7 millisecond pulsars with ISO. The pulsar 
PSR B1257+12 is orbited by three planets, and other millisecond pulsars may be orbited by 
dust disks that represent planets that failed to form or their residue. We do not detect any 
infrared emission from the 7 pulsars in our sample, and typical upper limits are 100 mJy. Using a 
simple model, we constrain the typical dust disk mass to be less than of order 100 M©, assuming 
that the heating of any putative dust disk would be coupled only weakly to the pulsar's emission. 
If the planets around PSR B1257+12 are composed largely of metals, our limits are probably 
an order of magnitude above plausible values for the disk mass in metals. Future observations 
with the Spitzer Space Telescope should be capable of probing into the range of disk masses that 
could plausibly give rise to planets. 

Subject headings: infrared: stars — planetary systems: protoplanetary disks — pulsars: general 



1. Introduction 

The first extrasolar planets discovered were 
found around the millisecond pulsar PSR B1257+12 
(Wolszczan & Frail 1992). The system consists 
of (at least) three planets, planet A with ap- 
proximately a lunar mass, planet B with M = 
4.3±0.2M e , and planet C with M = 3.9±0.2M ffi 
(Konacki & Wolszczan 2003). Although planetary 
systems around main-sequence stars had been long 
anticipated and numerous such systems have been 
found since, pulsar planetary systems were unex- 
pected. It was assumed that any planets orbiting 
the pulsar progenitor would have become gravita- 
tionally unbound in the supernova that produced 
the pulsar. 

Various mechanisms have been proposed for the 
formation of these planets (Phinney & Hansen 
1992; Miller & Hamilton 2001; Hansen 2002), but 
all generally rely on an accretion disk around the 
pulsar within which the planets form. Millisecond 
pulsars are a class of pulsars that, subsequent to 
their formation, undergo an episode of mass accre- 



tion from a companion (Wijnands & van der Klis 
1998). This process is thought to occur via an 
accretion disk, which transfers angular momen- 
tum to the pulsar as well, thereby spinning it up. 
Various mechanisms exist to shut down the ac- 
cretion (e.g., evolution of the companion), but, if 
the accretion is not 100% efficient, the millisec- 
ond pulsar will be left with an orbiting disk of 
material. Such a residual accretion disk is a nat- 
ural location for the formation of planets. Even 
if planets form, the formation process may leave 
a debris disk. Gozdziewski et al. (2003) investi- 
gated the long-term stability of a debris disk in 
the PSR B1257+12 system, finding a stable zone 
outside 1 AU. 

Since the discovery of planets around PSR 
B1257+12, a planet has also been found around 
the pulsar PSR B1620-26 (Backer, Foster, & Sail- 
men 1993; Thorsett, Arzoumanian, & Taylor 1993; 
Rasio 1994; Joshi & Rasio 1997; Thorsett et al. 
1999; Ford et al. 2000; Sigurdsson et al. 2003; 
Richer et al. 2003), in the globular cluster M4. In 
contrast to the planets orbiting PSR B1257+12, 
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which are thought to have formed in situ, the 
planet orbiting PSR B1620— 26 is thought to have 
been acquired during a dynamical exchange within 
the globular cluster. 

Pulsar planetary systems offer valuable in- 
sights, even if their total number is unlikely ever to 
approach the number of planetary systems around 
main-sequence stars. Taken together the two pul- 
sar planetary systems already indicate that plan- 
ets can form and exist in a wide variety of environ- 
ments. The presence of terrestrial mass planets 
around PSR B1257+12 suggests that terrestrial 
planets may be widespread, a hypothesis to be 
tested by future space missions such as Kepler 
and the Terrestrial Planet Finder (TPF). Plan- 
ets orbiting main-sequence stars near the Sun are 
found almost exclusively around stars with solar- 
or super-solar metallicities (Gonzalez 1997; San- 
tos et al. 2001), which has led to the belief that 
only stars with high metallicities can host planets. 
In contrast, the planet around PSR B1620— 26, 
if it was acquired during a dynamical exchange, 
probably has existed for a substantial fraction of 
the age of the globular cluster M4. This is a 
low-metallicity globular cluster, suggesting that 
planets can form in low-metallicity environments. 

Although the notion that planets can form in 
a residual accretion disk is plausible, no such ex- 
amples of residual accretion disks are known. The 
presence of a planet (or stellar companion) can 
be inferred using traditional pulsar timing tech- 
niques from the periodic advance and delay of the 
arrival time of the pulsar's pulse, due to the pul- 
sar's reflex motion. A relatively uniform disk of 
material would produce little reflex motion and 
therefore would remain undetected by these tra- 
ditional techniques. Detecting dust disks around 
millisecond pulsars not only would elucidate the 
late stages of millisecond pulsar "spin up" and 
planet formation, it would be a new probe of the 
local environments around millisecond pulsars. 

A modest number of unsuccessful searches 
for infrared emission from dust disks around 
millisecond pulsars have been conducted. Fig- 
ure 1 summarizes the current situation using 
PSR B1257+12 as an example. The limits for 
other pulsars are similar. 

This paper reports 60 and 90 /xm observations 
of 7 pulsars with the ISOPHOT instrument on- 
board the ISO satellite. In §2 we describe the ob- 
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Fig. 1. — A summary of current limits on infrared 
dust emission from around millisecond pulsars us- 
ing PSR B1257+12 as an example. Shown are 
limits placed by 10 /zm imaging by (Foster & Fis- 
cher 1996, FF96), IRAS limits determined by Fos- 
ter & Fischer (1996), limits from 800 /im imaging 
by (Phillips & Chandler 1994, PC94), and limits 
from 850 [im imaging by (Greaves & Holland 2000, 
GHOO). Also shown are the expected confusion 
limits for the future observations with the Spitzer 
Space Telescope at 70 and 160 /im. This pulsar 
was not observed as part of this ISO program. 
The solid curve shows the expected emission from 
a 300 M ffi dust disk composed of dust particles 
0.1 /im in size and heated by 1% of the spin-down 
luminosity of the pulsar while the dotted curve 
shows the expected emission from a 10 M s disk 
composed of 0.1 fim dust particles and heated by 
0.1% of the spin-down luminosity (Foster & Fis- 
cher 1996). Upper limits on infrared emission for 
other millisecond pulsars are similar. 
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servations and present our results and in §3 we 
describe how our results constrain the presence of 
dust disks around millisecond pulsars and present 
our conclusions. 

2. Observations and Data Analysis 

We compiled a list of millisecond pulsars known 
prior to 1994 August and with distances less than 
1 kpc. Distances are estimated from the Taylor 
& Cordes (1993) model and should be accurate 
to approximately 25%. Most of these millisecond 
pulsars lie at high Galactic latitudes. 

Of these, seven were observed with the ISOPHOT 
instrument (Lemke et al. 1996) onboard the ISO 
satellite (Kessler et al. 1996) between 1996 August 
and 1997 May. Table 1 summarizes the observing 
details; we also report the distance to each pulsar 
and, anticipating later discussion, whether or not 
it is a binary and its spin-down luminosity. 1 All 
of the observations used the P32 observing mode 
with the C100 detector. In this mode the space- 
craft was commanded to cover a series of raster 
pointings around the nominal pulsar position. At 
each raster pointing an internal chopper pointed 
the beam toward 13 adjacent sky positions. The 
throw of the chopper was larger than the offset 
between raster pointings. The result was that, 
in general, an individual sky position within the 
raster was observed multiple times or oversampled. 
Before and after each observation of a pulsar, an 
internal calibration source was observed. 

The analysis of the pulsar observations largely 
followed the standard ISOPHOT analysis pipeline. 
The key difference was the amount of "deglitch- 
ing" performed. Glitches result from cosmic rays 
striking the detector or secondary electrons pro- 
duced by spacecraft materials struck by primary 
cosmic rays. Failure to remove glitches can cor- 
rupt later calibration of all data, not just of the 
portion containing the glitches. The standard 
ISOPHOT analysis pipeline removes glitches but 
does so without making use of the redundancy im- 
plicit in the oversampled P32 observations. 

Deglitching proceeded in the following fash- 
ion. Within each spacecraft pointing the chopper 

1 The spin-down luminosity of a pulsar is a measure of its 
energy-loss due to magnetic dipole radiation and is given 
by L = IQQ, where / is its moment of inertia and Q is its 
rotation frequency. 



would sweep past a particular sky position mul- 
tiple times (typically 3-5 times). For each sky 
position, the median signal level was determined, 
then subtracted from all observations at that sky 
position. The observations from all sky positions 
were then combined to form a signal strength his- 
togram. A signal strength threshold was specified, 
and signals above this level were eliminated. Typ- 
ically 3%-10% of the signals were eliminated in 
this stage. Depending upon the number of chop- 
per sweeps per spacecraft pointing and deglitching 
prior to this stage, the median signal strength per 
sky position could not always be determined ac- 
curately. Thus, additional manual deglitching was 
done to remove any remaining outlier signals. Our 
use of the observations of the internal calibration 
sources followed the standard ISOPHOT analysis 
pipeline. 

After deglitching and calibration using the 
internal calibration sources, mapping was done 
within the ISOPHOT Interactive Analysis pack- 
age. Measurements from the individual detector 
pixels were co-added to form a sky image, with the 
contributions from the individual detector pixels 
weighted by their distances from the image pix- 
els. Doing so takes into account the beam profile 
falling on each detector pixel. We also employed a 
median flat field, which has the effect of reducing 
substantially our sensitivity to any extended emis- 
sion in the field. As we are attempting to detect 
point sources, we regard this reduced sensitivity 
to extended emission as unimportant. 

In no case have we identified a source at the 
location of a pulsar. Utilizing the inner quarter of 
the image, we determined the rms noise level. We 
take our upper limits to be 3 times this rms noise 
level. Table 2 summarizes the upper limits. 

3. Discussion and Conclusions 

We have not detected infrared emission associ- 
ated with any of the pulsars observed with ISO. 
Other infrared and sub-millimeter observations of 
millisecond pulsars have been conducted, and all 
of these have yielded only upper limits as well. 
Those observations most relevant to our sample of 
millisecond pulsars are those by Foster & Fischer 
(1996) at 10 fiia and Greaves & Holland (2000) 
at 850 /im. Foster & Fischer (1996) also utilized 
IRAS observations to obtain upper limits on the 
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Table 1 
Pulsars Observed 



Spin-down On-source 



Name 


Binary? 


Distance 
(kpc) 


Luminosity 

(Lo) 


Wavelength 
(//m) 


P32 Raster 


Time 

00 


PSR J0034-0534 


Y 


0.98 


10 


90 


3 


X 


8 


1402 


PSR J1640+2224 


Y 


1.19 


0.88 


60 


3 


X 


6 


1012 










90 


3 


X 


8 


848 


PSR J1730-2304 




0.51 


< 0.35 


60 


3 


X 


6 


1590 










90 


3 


X 


8 


1232 


PSR B1855+09 


Y 


0.91 


1.1 


60 


3 


X 


6 


1012 


PSR J2124-3358 




0.25 


1.7 


60 


3 


X 


6 


1012 










90 


3 


X 


8 


848 


PSR J2145-0750 


Y 


0.50 


< 0.048 


60 


3 


X 


6 


1590 










90 


3 


X 


8 


848 


PSR J2322+2057 




0.78 


0.62 


60 


3 


X 


6 


1804 










90 


3 


X 


8 


1402 



Table 2 

Dust Disk IR Emission Upper Limits 



Name Flux Density 

(mJy) 



60 ^.m 



PSR J1640+2224 


80 


PSR J1730-2304 


35 


PSR B1855+09 


190 


PSR J2124-3358 


100 


PSR J2145-0750 


100 


PSR J2322+2057 


59 


90 fim 


PSR J0034-0534 


48 


PSR J1640+2224 


59 


PSR J1730-2304 


140 


PSR J2124-3358 


55 


PSR J2145-0750 


73 


PSR J2322+2057 


39 



Note. — Upper limits are 3a. 
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infrared emission from their sample of pulsars. As 
Figure 1 shows, the upper limits set by IRAS are 
typically well above the limits set by our ISO ob- 
servations. Moreover, there is unfortunately lit- 
tle overlap between these three samples of pulsars 
(those whose observations are reported here, Fos- 
ter & Fischer 1996, and Greaves & Holland 2000). 
Most of the pulsars that have been observed be- 
tween 10 and 850 fim have been observed at only 
one or two wavelengths. 

Foster & Fischer (1996) developed a model for 
the infrared emission from a dust disk around a 
millisecond pulsar. Their model assumes that the 
disk consists of particles of a uniform radius a 
heated by a fraction / s( j of the pulsar's spin-down 
luminosity L s( j. The total mass of the disk is ma. 
While the model is simplistic — an actual dust disk 
presumably consists of particles with a range of 
sizes, the heating mechanism is left unspecified, 
non-equilibrium effects such as stochastic heating 
are ignored, and the impact of any stellar com- 
panions (see Table 1) on the disk are ignored — we 
believe that this simplicity is justified given the 
uncertainties of the heating mechanism and of the 
environs of a millisecond pulsar. 

In this model, for / s( j ~ 1%, typical dust tem- 
peratures are predicted to be T w 10-50 K for 
disks having nid ~ 100 M ffi and a ~ 1 /im and 
heated by a pulsar with L s( j ~ 1 L Q . These 
temperatures are similar to the lower temperature 
range used by Koch-Miramond et al. (2002) and 
considerably lower than those assumed (w 150 K) 
by Phillips & Chandler (1994), who estimated 
disk temperatures by scaling from observations 
of T Tauri stars. The lower temperatures result 
from our assumption of a weaker coupling between 
the pulsar's spin-down luminosity and the disk. 
Phillips & Chandler (1994) considered disk tem- 
perature to be a major uncertainty in convert- 
ing from measured flux densities to inferred disk 
masses. Accordingly, our assumption of a weaker 
coupling means that larger disk masses can be 
tolerated without violating the observational con- 
straints. 

Given the paucity of data, it is not possible, in 
general, to constrain all three parameters of this 
model with the existing observations. We there- 
fore adopt an approach in which we infer limits 
on two parameters of the Foster & Fischer (1996) 
model for fiducial values of the third parameter. 



Here, as an example, we consider the millisecond 
pulsar PSR J0034-0534 (Bailes et al. 1994) which 
has a probable white dwarf companion, is at a dis- 
tance of 1 kpc, and has a spin-down luminosity 
of 10 L Q . Greaves & Holland (2000) placed a 2a 
limit of 3.7 mJy at 850 /im, and we place a 2a 
limit of 50 mJy at 90 /jm. Figure 2 shows the 
allowed region of the disk mass-grain size plane 
given these observational limits and an assumed 
heating efficiency of / s d = 1%. 



o 
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Fig. 2. — Allowed disk masses and grain sizes for 
a dust disk orbiting PSR J0034-0534. The cross- 
hatched region indicates disk masses and grain 
sizes that do not violate the observational limits 
at 90 and 850 /im. A heating efficiency of / sc i = 1% 
has been assumed. 

Allowed regions in the m^-a plane occur for one 
of two possible reasons. First, the peak of the dust 
disk emission may appear shortward of 90 /im, 
where no constraints exist for this pulsar, with 
the Rayleigh-Jeans tail of the emission falling be- 
low the two measured values. This region is to 
the lower left in Figure 2. Second, the peak of the 
emission may appear between 90 fim and 850 /im, 
but with a magnitude comparable to that mea- 
sured at 90 fim so that the Rayleigh-Jeans tail 
again does not violate the 850 /im limit while 
the Wien tail of the emission does not violate the 
90 /im limit. This region is to the lower right in 
Figure 2. Obviously, a lower value of / s a would 
produce larger allowed regions in the nid-a plane. 
Larger allowed regions would also exist for other 
pulsars (Table 1) with smaller spin-down luminosi- 
ties. 

We conclude that, with the current observa- 
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tional constraints and the assumption of a fairly 
weak energy coupling between pulsars and disks, 
dust disks of order 100 M ffi easily could exist 
around millisecond pulsars. Koch-Miramond et 
al. (2002) have reached similar conclusions based 
on 15 and 90 fim observations. To the extent that 
such disks would be uniform, they would also es- 
cape detection from traditional pulsar timing tech- 
niques. Pulsar timing techniques utilize the ad- 
vance or delay of the pulse arrival time resulting 
from the pulsar's reflex motion to detect planetary 
or stellar companions. A relatively uniform dust 
disk would produce little reflex motion. 2 

The current limits on dust disk masses are far 
larger than the mass of the disk thought to have 
produced the planets around PSR B1257+12. The 
minimum combined mass of the two larger planets 
in that system is 8.2 M e (Konacki & Wolszczan 
2003). Assuming that these planets are composed 
largely of metals, we expect that any dust mass 
prior to the planets' formation would be compa- 
rable in magnitude. Indeed, Hansen (2002) has 
shown how an initial disk of mass 0.1-10~ 3 M Q 
containing of order 10 M e in metals could form 
a system similar to that orbiting PSR B1257+12. 
(See also Miller & Hamilton 2001 and Figure 1.) 
For those pulsars orbited by stellar companions, 
the companions will introduce regions of limited 
orbital stability within the disks, potentially im- 
plying even smaller expected disk masses. We 
conclude that current observational limits on dust 
disk masses are at least an order of magnitude 
above plausible values. 

The mid- to far-infrared detectors (24, 70, 
and 160 ^m) on the Spitzer Space Telescope 
should have sensitivities some 1-2 orders of mag- 
nitude better than the limits we report here. At a 
minimum, we expect that future observations with 
the Spitzer Space Telescope may require more so- 
phisticated modelling of disks, including the possi- 
ble effects of stellar companions. If future observa- 
tions with the Spitzer Space Telescope do not de- 
tect infrared emission from dust disks around mil- 
lisecond pulsars, the resulting mass limits should 
be in the range of 10 M 9 , sufficient to begin plac- 

2 One possible exception to this conclusion would be a dust 
disk illuminated directly by the pulsar's beam. In this case 
the relativistic particle flow from the pulsar beam poten- 
tially could ionize portions of the disk and induce plasma 
propagation delays which might be detectable. 



ing stringent constraints on their existence, or the 
temperatures of any pulsar dust disks must be no 
more than a few Kelvin. 
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